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ABSTRACT 


Conventional  separation  techniques  such  as  filters  cannot  be  used  in  a  scenario  where  a 
weak  signal  is  embedded  within  a  stronger  signal  in  the  same  frequency  band.  An  adap¬ 
tive  filter  approach  to  recover  a  weak  narrowband  communication  signal  embedded  in  a 
strong  wideband  broadcast  signal  and  additive  noise  is  investigated  by  taking  advantage  of 
the  known  characteristics  of  the  broadcast  signal  standard.  The  weak  signal  considered  is 
a  four-quadrature  amplitude  modulation  (QAM)  communication  signal,  and  the  stronger 
signal  is  a  Digital  Terrestrial  Multimedia  Broadcasting  (DTMB)  signal,  which  is  the  digital 
television  standard  for  People’s  Republic  of  China.  The  results  show  that,  while  the  extrac¬ 
tion  of  the  weak  signal  is  possible,  the  characteristics  of  the  transmission  channel  must  be 
very  accurately  estimated  for  the  scheme  to  achieve  a  reasonable  error  rate. 


V 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


VI 


Table  of  Contents 


1  Introduction  1 

1.1  Motivation .  1 

1.2  Past  Work .  1 

1.3  Approach .  2 

1.4  Thesis  Organization .  3 

2  Theoretical  Background  5 

2. 1  Adaptive  Noise  Cancellation .  5 

2.2  Structure  of  DTMB  Signal .  6 

2.3  Channel  Estimation  for  DTMB  Signal .  6 

3  Methodology  9 

3.1  Schematic  of  Overall  Implementation .  9 

3.2  Preprocessing  for  Adaptive  Noise  Cancellation .  9 

3.3  Adaptive  Noise  Cancellation .  10 

3.4  Simulink  Implementation .  11 

3.5  Simulations .  12 

4  Results  and  Discussions  15 

4.1  Need  for  Multiple  Trials .  15 

4.2  Effect  of  Number  of  Coefficients  in  Channel  Estimate .  16 

4.3  Effect  of  Error  in  Channel  Coefficients .  17 

5  Conclusions  25 

5.1  Summary .  25 

5.2  Recommendations  for  Euture  Work .  25 


Appendix  A  Simulink  Models 


27 


Appendix  B  MATLAB  Code  31 

B.l  Code  to  Run  Simulation .  31 

B.2  Code  for  Initialization .  34 

B.3  Code  to  Enable  Parallel  Computing .  37 

Appendix  C  Supplementary  Plots  43 

List  of  References  45 

Initial  Distribution  List  49 


viii 


List  of  Figures 


Figure  1.1  Plot  showing  a  weak  eo-ehannel  eommunieation  signal  embedded 

within  a  strong  television  broadeast  signal .  2 

Figure  1.2  Areas  of  the  world  where  ATSC,  DVB-T  and  DTMB  are  used.  .  3 

Figure  2. 1  Typieal  strueture  of  a  basie  adaptive  noise  eaneeler.  .  5 

Figure  2.2  Strueture  of  the  DTMB  signal  frame .  6 

Figure  2.3  Frame  headers  of  various  length .  7 

Figure  3.1  Sehematie  of  overall  setup .  9 

Figure  3.2  Preproeessor  for  adaptive  noise  eaneellation .  10 

Figure  3.3  Adaptive  noise  eaneeler  bloek .  11 

Figure  3.4  Simulink  implementation  of  the  preproeessor.  .  12 

Figure  4. 1  Plot  of  SER  against  INR  with  42  000  symbols  transmitted  per  INR.  17 

Figure  4.2  Plot  of  SER  against  INR  with  42000  symbols  transmitted  per  INR 

over  100  trials .  19 

Eigure  4.3  Plot  of  mean  SER  against  INR  with  different  number  of  eoeffieients 

k  in  the  estimated  CIR  averaged  over  100  trials .  20 

Eigure  4.4  Plot  of  mean  SER  against  INR  for  different  number  of  eoeffieients  k 

in  the  estimated  CIR  averaged  over  100  trials  with  error  bars.  .  .  21 

Eigure  4.5  Typieal  plot  of  estimated  CIR  fi(n) .  22 

Eigure  4.6  Plot  of  mean  SER  against  INR  for  different  errors  in  estimated  ehan- 

nel  eoeffieients  averaged  over  100  trials .  22 

Eigure  4.7  Plot  of  mean  SER  against  INR  for  different  errors  in  estimated  ehan- 

nel  eoeffieient  values  averaged  over  100  trials  with  error  bars.  .  .  23 


IX 


Figure  A.  1  Overall  Simulink  model .  27 

Figure  A.2  Simulink  model  for  eommunieation  signal  modulator .  28 

Figure  A. 3  Simulink  model  for  DTMB  modulator. .  28 

Figure  A. 4  Simulink  model  for  transmission  ehannel .  28 

Figure  A. 5  Simulink  model  for  preprocessor  and  adaptive  noise  canceler.  .  .  29 

Figure  A. 6  Simulink  model  for  the  DTMB  demodulator  and  the  channel  esti¬ 
mator  within  the  DTMB  demodulator. .  30 

Figure  A. 7  Simulink  model  for  the  first  communication  signal  demodulator.  .  30 

Figure  A. 8  Simulink  model  for  the  second  communication  signal  demodulator.  30 

Figure  C.l  Plot  showing  distribution  of  the  100  observations  of  SER  for  each 

INR  value  for  different  number  of  coefficients  k  in  h{n) .  43 

Figure  C.2  Plot  showing  distribution  of  the  100  observations  of  SER  for  each 
INR  value  for  different  percentages  of  error  in  the  channel  coeffi¬ 
cients .  44 


X 


List  of  Tables 


Table  3.1  Simulation  parameters  for  television  broadcast  signal .  13 

Table  3.2  Simulation  parameters  for  the  communication  signal .  13 

Table  3.3  Characteristics  for  the  multipath  channel .  14 

Table  3.4  Characteristics  for  the  NLMS  adaptive  filter. .  14 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


List  of  Acronyms  and  Abbreviations 


ATSC 

AWGN 

CIR 

DTMB 

DVB-T 

FFT 

INR 

NLMS 

OFDM 

SER 

SIR 

SNR 

QAM 


Advanced  Television  Systems  Committee 
additive  white  Gaussian  noise 
channel  impulse  response 
Digital  Terrestrial  Multimedia  Broadcasting 
Digital  Video  Broadcasting-Terrestrial 
fast  Fourier  transform 
interference-to-noise  ratio 
normalized-least-mean-square 
orthogonal  frequency-division  multiplexing 
symbol  error  ratio 
signal-to-interference  ratio 
signal-to-noise  ratio 
quadrature  amplitude  modulation 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


XIV 


Executive  Summary 


Consider  a  scenario  where  a  weak  signal  such  as  a  communication  signal  is  embedded 
within  a  strong  signal  such  as  a  television  broadcast  signal.  Both  signals  overlap  the  same 
frequency  band  and  cannot  be  separated  by  filtering  in  the  frequency  domain.  We  inves¬ 
tigate  an  adaptive  filter  approach  to  recover  the  weak  narrowband  communication  signal 
embedded  in  the  strong  wideband  television  broadcast  signal  and  additive  white  Gaussian 
noise  (AWGN)  by  taking  advantage  of  the  known  characteristics  of  the  television  broadcast 
standard.  If  successful,  the  technique  could  be  developed  to  allow  the  military  to  create  a 
covert  communication  channel  or  offer  protection  against  jamming.  In  civilian  applica¬ 
tions,  this  could  allow  reuse  of  the  frequency  band  now  used  by  television  broadcasters, 
creating  more  bandwidth  available  for  transmission. 

Cristi,  Fargues,  and  Hagstette  have  previously  studied  a  similar  problem  in  the  context  of 
a  weak  co-channel  communication  signal  embedded  within  a  strong  terrestrial  television 
broadcast  signal  [1].  They  identified  three  possible  schemes  to  separate  the  signals,  namely 
blind  source  separation,  using  multiple  antennas  and  Kalman  filtering,  and  using  adaptive 
noise  cancellation.  The  first  two  were  studied  by  Hagstette  [2],  [3]  and  Attique  [4].  The 
adaptive  noise  cancellation  approach  is  the  focus  of  this  study. 

In  these  past  studies,  the  television  broadcast  signal  types  considered  followed  the  Advanced 
Television  Systems  Committee  (ATSC)  standard,  used  in  North  America,  and  the  Digital 
Video  Broadcasting-Terrestrial  (DVB-T)  standard,  used  in  Europe,  Africa,  India,  and  oth¬ 
ers.  Recently,  Lai  developed  a  Simulink  model  to  generate  signals  following  the  Digital 
Terrestrial  Multimedia  Broadcasting  (DTMB)  standard  used  in  People’s  Republic  of  China 
[5].  Being  a  more  recent  standard,  DTMB  uses  more  advanced  error  correction  and  is  more 
robust  to  Doppler  shift  distortion  than  previously  adopted  standards  [6],  [7].  In  this  study, 
the  signals  considered  are  a  four-quadrature  amplitude  modulation  (QAM)  communication 
signal  and  a  DTMB  television  broadcast  signal.  We  assume  the  communication  signal  and 
DTMB  signal  transmitters  are  not  co-located.  The  DTMB  signal  transmission  channel  has 
multiple  paths,  while  the  communication  signal  transmission  channel  has  a  single  path  and 
AWGN. 


XV 


We  implement  a  Simulink  model  of  the  proposed  approach  and  present  the  schematic  of  the 
overall  setup  in  Figure  1 .  The  details  of  the  preprocessor  and  adaptive  noise  canceler  are 
shown  in  Figures  2  and  3.  The  DTMB  television  broadcast  signal  has  well-known  charac¬ 
teristics,  good  error  correction  and  strong  power.  Known  standard  characteristics  are  taken 
advantage  of  to  extract  the  weak  communication  signal  using  an  adaptive  noise  cancella¬ 
tion  scenario.  The  proposed  signal  extraction  process  has  two  phases.  First,  a  preprocessor 
reconstructs  the  DTMB  signal  and  minimizes  the  contribution  of  the  DTMB  signal  in  the 
received  combined  noisy  signal.  Second,  an  adaptive  noise  canceler  is  designed  to  recover 
the  communication  signal  using  the  reconstructed  DTMB  signal  as  a  reference. 

Results  show  that  the  adaptive  noise  canceler  requires  a  very  accurate  channel  estimation 
process  to  recover  the  weak  communication  signal.  Results  also  show  that  symbol  error 
ratio  (SER)  levels  below  10^^  may  be  obtained  with  up  to  20%  error  in  estimated  trans¬ 
mission  channel  coefficients  values  when  an  exact  length  of  the  transmission  channel  is 
selected. 
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CHAPTER  1: 
Introduction 


The  motivation  for  the  thesis,  past  work,  approach  followed,  and  thesis  organization  are 
discussed  in  this  chapter. 

1.1  Motivation 

In  this  thesis,  we  consider  a  scenario  where  a  weak  signal  such  as  a  communication  sig¬ 
nal  is  embedded  within  a  much  stronger  signal  such  as  a  television  broadcast  signal  in  the 
same  frequency  band,  as  shown  in  Figure  1.1.  Since  both  signals  overlap  the  same  fre¬ 
quency  band,  they  cannot  be  separated  by  filtering  in  the  frequency  domain.  An  adaptive 
filter  approach  to  recover  a  weak  narrowband  communication  signal  embedded  in  a  strong 
wideband  television  broadcast  signal  and  additive  noise  is  investigated  by  taking  advantage 
of  the  known  characteristics  of  the  television  broadcast  signal  standard.  This  could  have 
both  military  and  civilian  applications  if  successful.  For  instance,  in  the  military,  this  could 
be  used  to  create  a  covert  communication  channel  or  offer  protection  against  jamming. 
In  civilian  applications,  this  could  allow  for  reuse  of  the  same  frequency  band  and  more 
bandwidth  available  for  transmission. 


1.2  Past  Work 

Radio  frequency  is  a  limited  resource.  To  maximize  the  use  of  this  limited  resource,  one  of 
the  techniques  used  by  communication  systems  is  frequency  reuse.  Often,  frequency  reuse 
results  in  co-channel  interference,  limiting  performance  [1].  There  has  been  much  research 
conducted  to  resolve  co-channel  interference  in  areas  such  as  cellular  communications  and 
satellite  communications  [2]-[6]. 

Cristi,  Fargues,  and  Hagstette  have  previously  studied  a  similar  problem  in  the  context  of 
a  weak  co-channel  communication  signal  embedded  within  a  strong  terrestrial  television 
broadcast  signal  [7].  They  identified  three  possible  schemes  to  separate  the  signals,  namely 
blind  source  separation,  using  multiple  antennas  and  Kalman  filtering,  and  using  adaptive 
noise  cancellation.  The  first  two  were  studied  by  Hagstette  [8],  [9]  and  Attique  [10].  The 
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Figure  1.1.  Plot  showing  a  weak  co-channel  communication  signal  (green)  embedded  within  a 
strong  television  broadcast  signal  (red). 

third  approach  of  using  adaptive  noise  cancellation  is  the  focus  of  this  study. 

In  these  previous  studies,  the  television  broadcast  signal  types  considered  followed  the 
Advanced  Television  Systems  Committee  (ATSC)  standard,  used  in  North  America,  and 
the  Digital  Video  Broadcasting-Terrestrial  (DVB-T)  standard,  used  in  Europe,  Africa,  In¬ 
dia,  and  others,  as  seen  in  Figure  1.2.  These  television  standards  were  developed  in  the 
1990s.  Recently,  Lai  developed  a  Simulink  model  to  generate  signals  following  the  Digital 
Terrestrial  Multimedia  Broadcasting  (DTMB)  standard  used  in  China  [11].  The  DTMB 
standard  was  adopted  as  a  national  standard  of  People’s  Republic  of  China  in  2007  [12]. 
Being  a  more  recent  standard,  DTMB  uses  more  advanced  error  correction  and  has  a  bet¬ 
ter  performance  with  respect  to  Doppler  shift  than  previously  adopted  standards  [13],  [14]. 
These  advantages  allow  for  easier  signal  reconstruction  during  demodulation  and  also  make 
the  adaptive  noise  cancellation  step  easier  to  implement. 

1.3  Approach 

In  this  study,  we  consider  a  four-quadrature  amplitude  modulation  (QAM)  signal  to  repre¬ 
sent  the  weak  communication  signal  and  a  DTMB  signal  to  represent  the  strong  signal.  A 
4— QAM  signal  was  selected  to  represent  a  typical  digital  communication  signal.  As  men- 
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Figure  1.2.  Areas  of  the  world  where  ATSC,  DVB-T  and  DTMB  are  used,  from  [15]. 

tioned  earlier,  the  DTMB  signal  is  chosen  because  its  model  is  readily  available  from  past 
work,  and  the  standard  is  known  to  have  better  performance  than  older  standards  [14].  We 
test  the  setup  with  a  simulation  model  in  Simulink. 

1.4  Thesis  Organization 

In  Chapter  2,  we  provide  the  theoretical  background  needed  to  understand  the  method  used 
in  this  thesis.  In  Chapter  3,  we  illustrate  the  method  and  its  implementation  in  Simulink. 
In  Chapter  4,  we  present  the  results  obtained.  In  Chapter  5,  we  provide  conclusions  and 
recommendations  for  future  work.  In  Appendixes  A  and  B,  we  show  the  MATLAB  code 
and  the  Simulink  models  used  in  the  simulation.  In  Appendix  C,  we  present  additional 
plots  to  supplement  the  information  provided  in  Chapter  4. 


3 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


4 


CHAPTER  2: 
Theoretical  Background 


In  this  chapter,  we  first  present  a  brief  introduction  to  adaptive  noise  cancellation.  Next,  we 
discuss  the  structure  of  a  DTMB  signal  and  channel  estimation  for  a  DTMB  signal. 


2.1  Adaptive  Noise  Cancellation 

Adaptive  noise  cancellation  has  been  used  to  remove  unwanted  interference  from  signals 
in  numerous  applications.  Some  of  the  earlier  applications  are  used  in  the  removal  of  un¬ 
wanted  interference  in  electrocardiograms,  canceling  noise  in  speech  signals,  or  removing 
antenna  sidelobe  interference.  The  unwanted  interference  could  be  the  60  Hz  from  the 
power  line  in  a  patient’s  electrocardiogram,  the  mother’s  own  electrocardiogram  in  a  fe¬ 
tus’s  electrocardiogram,  engine  noise  in  a  pilot’s  speech,  or  an  unwanted  signal  coming 
from  the  sidelobes  in  an  antenna  [16],  [17]. 

A  typical  structure  of  a  basic  adaptive  noise  canceler  is  shown  in  Figure  2.1.  The  sensor 
input  has  the  signal  s{n)  and  an  uncorrelated  noise  vi{n)  corrupting  the  signal.  We  desire  to 
remove  the  noise  v\ (n)  from  the  signal.  A  second  sensor  provides  a  reference  for  the  noise 
V2{n),  which  is  correlated  with  the  noise  vi{n).  The  aim  of  the  filter  is  to  take  advantage  of 
the  correlation  between  vi{n)  and  V2{n)  to  produce  an  output  y{n)  such  that  y{n)  =  vi{n). 
After  subtraction,  the  error  signal  e{n)  will  be  the  signal  s{n)  [16]. 


Figure  2.1.  Typical  structure  of  a  basic  adaptive  noise  canceler,  after  [18]. 
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2.2  Structure  of  DTMB  Signal 

A  DTMB  signal  is  transmitted  in  frames  and  is  based  on  orthogonal  frequency-division 
multiplexing  (OFDM).  The  most  fundamental  frame  is  the  signal  frame.  As  shown  in 
Figure  2.2,  the  signal  frame  is  made  up  of  a  frame  header  and  a  frame  body  with  a  baseband 
symbol  rate  of  7.56  mega  symbols  per  second.  The  frame  body  has  exactly  3780  symbols. 
The  frame  header,  however,  could  have  a  length  of  420,  595,  or  945  symbols,  depending  on 
applications.  In  general,  longer  frame  headers  have  a  lower  throughput  but  higher  resilience 
to  adverse  conditions.  In  Figure  2.3,  we  can  see  the  structure  for  the  frame  header  of 
different  lengths.  Frame  headers  with  a  length  of  420  and  945  symbols  have  a  pre-amble, 
a  pseudo-noise  sequence,  and  a  post-amble.  Frame  headers  with  a  length  of  595  symbols 
have  only  pseudo-noise  sequence  and  no  pre-  and  post- ambles.  In  particular,  note  that  the 
pseudo-noise  sequences  can  be  used  as  pilot  symbols  to  perform  channel  equalization  as 
they  are  known  at  the  receiver  side  [13],  [19]. 


2.3  Channel  Estimation  for  DTMB  Signal 

In  this  section,  we  focus  on  one  aspect  of  DTMB  modulation  and  demodulation — channel 
estimation.  Transmission  channel  characteristics  have  a  direct  impact  on  signals  propagat¬ 
ing  through  them.  We  can  model  the  channel  impact  as  a  linear  filter  with  impulse  response 
h{n).  The  impulse  response  h{n)  is  known  as  the  channel  impulse  response  (CIR).  We  want 
to  undo  the  impact  of  the  channel  on  the  signal  at  the  receiver  so  that  we  can  correctly  inter¬ 
pret  the  signal.  The  removal  of  the  transmission  channel  distortion  on  the  signal  is  known 
as  channel  equalization.  Typically,  the  first  step  to  channel  equalization  is  to  estimate  the 
CIR,  that  is,  channel  estimation. 

A  DTMB  signal  is  transmitted  in  frames,  as  mentioned  in  Section  2.2.  There  are  pseudo¬ 
noise  sequences  within  the  frame  header  that  are  known  to  the  receiver.  The  receiver  can 
use  these  predetermined  symbols  as  pilot  symbols  to  perform  channel  estimation  [20],  [21]. 


Frame  header 


Frame  body 


Figure  2.2.  Structure  of  a  DTMB  signal  frame,  from  [13]. 
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Pre-amble  (82) 


Pseudo-noise  sequenee  (255) 


Post-amble  (83) 


(a)  Frame  header  of  length  420  symbols. 


Pseudo-noise  sequenee  (595) 

(b)  Frame  header  of  length  595  symbols. 


Pre-amble  (217) 


Pseudo-noise  sequenee  (511) 


Post-amble  (217) 


(c)  Frame  header  of  length  945  symbols. 

Figure  2.3.  Frame  headers  of  various  length.  The  number  in  parenthesis  denotes  the  number  of 
symbols,  after  [13],  [19]. 


The  estimated  channel  frequency  response  H{f)  can  be  calculated  by 


YpiUf) 

Xpilotif) 


(2.1) 


where  XpHotif)  is  the  transmitted  pilot  sequence  in  the  frequency  domain  and  YpHatif)  is 
the  received  pilot  sequence  in  the  frequency  domain. 

The  estimated  CIR  h{n)  can  then  be  found  from 


h{n)  =  IFFT  [H{f) 


(2.2) 


where  IFFT[.]  denotes  the  inverse  fast  Fourier  transform  (FFT)  operation. 

The  calculation  of  H{f)  using  (2.1),  however,  runs  a  risk  of  division  by  zero,  particularly 
when  the  quantity  XpHotif)  is  small. 


To  improve  the  numerical  stability,  we  note  that 

1  1 


V* 

pilot 


if) 


Xpiiotif)  Xpaorif)x;pjf) 

^UAf) 

\Xpilotif)\^ 


(2.3) 


V* 

pilot 


if) 


\Xpiiotif)\^  +  S 


1 


where  is  the  eomplex  eonjugate  of  XpUotif)  and  e  is  a  small  number. 

The  small  number  e  added  to  the  denominator  of  (2.3)  is  designed  to  improve  the  numerieal 
stability  of  the  division. 


Substituting  (2.3)  into  (2.1),  we  get 


\Xpiiotif)\^  +  £ 


(2.4) 


By  using  (2.4)  to  ealeulate  H{f),  we  get  improved  numerieal  stability  in  the  estimation  of 
the  transmission  ehannel  frequeney  response,  espeeially  when  the  quantity  XpUotif)  is  very 
small.  This  in  turn  ensures  the  estimated  CIR  is  numerieally  better  approximated. 
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CHAPTER  3: 
Methodology 


The  adaptive  noise  cancellation  process  implemented  to  extract  the  communication  signal 
embedded  in  the  television  broadcast  signal  is  discussed  in  this  chapter.  The  methodology 
was  developed  by  Cristi  and  used  in  this  thesis  with  permission  [7]. 

3.1  Schematic  of  Overall  Implementation 

In  Figure  3.1,  we  show  a  schematic  of  the  overall  setup.  On  the  left,  we  generate  a  4-QAM 
communication  signal  and  a  DTMB  television  broadcast  signal.  The  communication  chan¬ 
nel  has  additive  white  Gaussian  noise  (AWGN)  and  the  television  broadcast  signal  passes 
through  a  transmission  channel.  The  received  signal  undergoes  preprocessing,  as  discussed 
in  Section  3.2,  followed  by  adaptive  noise  cancellation,  as  discussed  in  Section  3.3.  We  then 
recover  the  communication  signal  and  compare  against  the  transmitted  version  to  compute 
the  symbol  error  ratio  (SER). 


3.2  Preprocessing  for  Adaptive  Noise  Cancellation 

In  Section  2.1,  we  noted  that,  for  the  basic  adaptive  noise  cancellation,  we  need  the  signal 
that  has  been  corrupted  with  noise  s{n)  -\-  v\ (n)  and  a  reference  noise  V2{n) .  In  our  context, 
the  television  broadcast  signal  is  considered  as  the  noise,  and  the  communication  signal  is 
the  signal-of-interest.  In  our  application,  we  only  have  the  combined  signal  where  the  com¬ 
munication  signal  is  embedded  within  the  television  signal.  We  now  discuss  the  process  to 


Comms 

signal 

(base¬ 

band) 


Video 

stream 


Figure  3.1.  Schematic  of  overall  setup. 
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obtain  a  reference  noise  signal  from  the  combined  signal. 

In  Figure  3.2,  we  show  the  preprocessor  for  adaptive  noise  cancellation.  The  received 
signal  contains  the  communication  signal  embedded  in  the  television  broadcast  signal.  In 
the  upper  path,  we  perform  a  first  cancellation  through  the  upper  path  using  the  estimated 
channel.  The  delay  in  the  upper  path  is  added  to  synchronize  the  signal  due  to  computa¬ 
tional  delays  in  the  lower  path.  In  the  lower  path,  the  received  signal  passes  through  the 
TV  demodulator  to  obtain  the  baseband  TV  signal  and  the  estimated  CIR.  The  baseband 
TV  signal  is  then  re-modulated  into  the  television  broadcast  signal.  The  reconstructed  tele¬ 
vision  broadcast  signal  is  filtered  through  the  estimated  channel  and  is  subtracted  from  the 
received  signal  (first  cancellation)  to  compute  the  signal  zi  (n) .  The  adaptive  noise  can¬ 
cellation  process,  therefore,  only  needs  to  cancel  the  remaining  noise  due  to  the  channel 
estimation  error.  The  reconstructed  television  broadcast  signal  is  also  passed  to  the  adaptive 
noise  canceler  as  the  reference  signal  Z2(«)- 


3.3  Adaptive  Noise  Cancellation 

In  Figure  3.3,  we  show  the  details  of  the  adaptive  noise  canceler.  Since  the  final  goal 
is  to  reconstruct  the  communication  signal  in  baseband,  the  two  signals  zi{n)  and  Z2(«) 
passed  from  the  preprocessor  are  first  demodulated  to  baseband.  Then  they  are  used  as  the 
input  and  reference  to  the  adaptive  noise  canceler.  At  the  output,  we  recover  the  baseband 
communication  signal. 


Received 

signal 


Input  of 
adaptive 
noise 
canceler 
zi(n) 


Reference 

signal 

ziin) 


Figure  3.2.  Preprocessor  for  adaptive  noise  cancellation,  after  [7]. 
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Recovered 

comms 

signal 


3.4  Simulink  Implementation 

We  implement  the  discussed  model  in  Simulink  with  the  differences  described  in  the  fol¬ 
lowing  sections.  The  detailed  Simulink  model  and  MATLAB  code  are  shown  in  Ap¬ 
pendixes  A  and  B. 

3.4.1  Replacement  of  Modulator  and  Demodulator 

The  reference  signal  Z2(^)  is  generated  by  demodulating  the  received  television  broadcast 
signal  and  remodulating  it  in  the  preprocessor  for  the  adaptive  noise  cancellation,  as  shown 
in  Figure  3.2.  We  can  still  recover  the  television  broadcast  signal  despite  the  interference 
caused  by  the  communication  signal  by  taking  advantage  of  the  error  correction  in  the 
television  broadcast  signal.  In  the  Simulink  implementation,  we  do  not  have  the  television 
broadcast  signal  demodulator  and  modulator  blocks  because  they  add  an  additional  level 
of  complexity  to  the  model.  By  not  implementing  the  demodulator  and  modulator,  we 
eliminate  potential  sources  of  error  coming  from  the  demodulator  and  modulator.  In  this 
work,  we  assume  the  television  broadcast  signal  can  be  recovered  with  100%  accuracy  due 
to  the  superior  error  correction  coding  in  DTMB  [14]. 

In  place  of  the  television  demodulator  and  modulator,  we  implement  only  the  channel  es¬ 
timator  of  the  demodulator.  The  actual  television  broadcast  signal  from  the  transmitter  is 
routed  into  the  preprocessor  to  replace  the  output  from  the  modulator.  See  Figure  3.4  for 
the  Simulink  implementation  of  the  preprocessor. 

3.4.2  Inclusion  of  Blocks  to  Calculate  Error  Rate 

For  the  purpose  of  quantifying  the  performance  of  the  scheme,  we  examine  the  SER  for  a 
given  signal-to-noise  ratio  (SNR).  We  include  in  the  Simulink  implementation  the  means 
to  calculate  the  SER  for  the  simulation  runs. 
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Input  of 
adaptive 
noise 
canceler 
zi(n) 


Channel 
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TV  signal 

from  - 

transmitter 

Figure  3.4.  Simulink  implementation  of  the  preprocessor. 
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3.5  Simulations 

In  this  study,  we  investigate  the  ability  to  reeover  the  low-power  eommunication  signal 
using  adaptive  noise  eancellation  under  the  following  conditions. 

1.  The  television  broadcast  signal  is  assumed  to  be  a  DTMB  signal.  The  weak  commu¬ 
nication  signal  is  a  4-QAM  signal. 

2.  The  television  broadcast  signal  and  communication  signal  are  random  binary  bit- 
streams. 

3.  The  television  broadcast  transmitter  and  the  communication  transmitter  are  not  co¬ 
located.  The  television  broadcast  signal  and  the  communication  signal  pass  through 
different  transmission  channels  to  arrive  at  the  receiver.  The  transmission  channel  of 
the  television  broadcast  signal  has  multipaths,  while  the  communication  signal  goes 
through  a  single  path  with  AWGN. 

4.  We  investigate  the  capabilities  of  the  proposed  approach  to  recover  the  weak  commu¬ 
nication  signal  as  the  AWGN  power  varies.  We  define  two  power  ratio  quantities — 
interference-to-noise  ratio  (INR)  and  signal-to-interference  ratio  (SIR) — given  by 

Weak  communication  signal  power 

INR  = -  (3.1) 

AWGN  power 
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and 


Television  broadcast  signal  power 
Weak  communication  signal  power 

Throughout  the  study,  we  maintain  a  fixed  SIR  level  equal  to  30  dB.  This  specific 
level  is  selected  to  be  consistent  with  past  studies  [8]. 


3.5.1  Simulation  Parameters  for  Television  Broadcast  Signal 

The  television  broadcast  signal  is  based  on  the  DTMB  standards  and  is  generated  using  the 
model  developed  by  Lai  [11].  The  simulation  parameters  selected  for  the  DTMB  signal  are 
shown  in  Table  3.1.  The  DTMB  signal  has  a  frame  header  of  420  symbols.  Within  the  frame 
header,  the  pseudo-noise  sequence  has  a  length  of  255  symbols  which  are  used  for  channel 
estimation  in  this  study.  The  total  frame  length  is  equal  to  4200  symbols.  The  sampling 
rate  is  equal  to  7.56  MHz  to  be  in  line  with  standard  broadcasting  applications  [11]. 

Table  3.1.  Simulation  parameters  for  television  broadcast  signal,  after  [11]. 


Description 

Value 

Sampling  rate  Fs 

7.56  MHz 

Length  of  frame 

4200  symbols 

Length  of  frame  header 

420  symbols 

Length  of  pilot  symbols 

255  symbols 

3.5.2  Simulation  Parameters  for  Communication  Signal 

The  communication  signal  is  a  4-QAM  communication  signal.  In  Table  3.2,  we  list  the 
simulation  parameters.  The  carrier  frequency  and  data  rate  of  the  signal  are  chosen  such 
that  they  are  an  integer  division  of  the  sampling  rate  Fs.  This  selection  was  made  so  that  the 
sampling  period  can  become  the  step  size  used  in  the  simulation.  Note  that  the  fundamental 
sample  time  would  be  reduced,  slowing  down  the  simulation,  if  the  carrier  frequency  and 
data  rate  are  not  integer  divisions  of  the  sampling  rate  Fs  [22] . 

Table  3.2.  Simulation  parameters  for  the  communication  signal. 


Description 

Value 

Carrier  frequency 

F 

10 

Data  rate 

F 

100 
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3.5.3  Simulation  Parameters  for  Multipath  Channel 

The  television  broadcast  signal  passes  through  a  multipath  channel  to  get  to  the  receiver.  To 
be  consistent  with  past  studies  [7],  the  multipath  channel  has  three  paths  with  characteristics 
shown  in  Table  3.3.  As  a  result  of  the  sampling  rate  and  the  length  of  the  longest  path, 
we  need  a  filter  length  of  17  to  model  the  multipath  transmission  channel. 


Table  3.3.  Characteristics  for  the  multipath  channel. 


Delay  (samples  /  ps) 

Power  (dB) 

Path  1 

0/0 

0 

Path  2 

10/1.32 

-10 

Path  3 

15/1.98 

-15 

3.5.4  Simulation  Parameters  for  Adaptive  Filter 

In  Table  3.4,  we  show  the  characteristics  of  the  adaptive  filter  used  in  the  adaptive  noise 
cancellation.  We  choose  a  normalized-least-mean-square  (NLMS)  adaptive  filter  for  its 
simplicity,  robustness,  and  ability  to  track  changes  in  the  signal  due  to  possible  changes 
in  the  transmission  channel  [23].  The  length  of  the  filter  is  selected  to  be  large  enough 
to  model  the  characteristics  of  the  multipath  channel.  The  step  size  is  selected  through 
trial  and  error  to  get  reasonable  performance  from  the  adaptive  noise  canceler.  The  filter 
weights  are  arbitrarily  initialized  to  zero  at  the  beginning  of  the  simulation. 

Table  3.4.  Characteristics  for  the  NLMS  adaptive  filter. 


Description 

Value 

Filter  length 

50 

Step  size 

0.0001 

Initial  filter  weight 

0 
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CHAPTER  4: 
Results  and  Discussions 


In  this  chapter,  we  present  the  simulation  results  obtained.  The  ability  of  the  entire  scheme 
to  recover  the  weak  communication  signal  is  evaluated  by  examining  the  SERs  across  dif¬ 
ferent  INRs. 


4.1  Need  for  Multiple  Trials 

Typically,  the  SER  is  expected  to  decrease  monotonically  as  the  signal-to-noise  ratio  in¬ 
creases  [24].  In  Eigure  4.1,  we  show  the  results  from  a  single  simulation  run  obtained  for 
a  transmission  length  of  42000  symbols  for  each  INR.  The  SER  does  not  decrease  mono¬ 
tonically  with  increasing  INR  as  expected.  Hence,  multiple  runs  are  needed  to  provide 
statistically  meaningful  results.  In  Eigure  4.2,  we  present  the  results  obtained  by  averaging 
over  100  trials.  In  Eigure  4.2a,  the  lower  and  upper  ends  of  the  error  bars  represent  the  5th 
and  95th  percentile  levels,  respectively.  We  note  that,  other  than  the  case  where  the  INR  is 
equal  to  5  dB,  there  are  more  than  50%  of  the  SER  values  that  are  equal  to  zero  for  each 
INR,  as  seen  in  Eigure  4.2b.  To  allow  us  to  plot  such  values  onto  the  logarithmic  axis,  they 
are  represented  as  10^^  in  the  plot.  The  non-zero  SER  values,  therefore,  skew  the  mean 
SER  upwards,  since  more  than  50%  of  the  SER  values  are  equal  to  zero.  We  note  that, 
since  more  than  50%  of  the  SER  values  are  equal  to  zero,  the  median  SER  is  zero. 

Increasing  the  number  of  trials  helped  to  produce  more  statistically  meaningfully  results.  It 
also,  however,  lengthened  the  time  required  to  produce  results  significantly.  Eor  instance, 
on  an  Intel  Xeon  computer  available  in  our  laboratory,  a  single  run  where  42000  symbols 
were  transmitted  took  about  0.25  h.  Generating  the  complete  SER  curve  over  a  range  of 
INRs  between  —  5  dB  and  50  dB  at  5  dB  interval  took  about  3  h. 

To  reduce  the  computational  impact,  we  used  the  Parallel  Computing  Toolbox  in  MATEAB 
to  enable  parallel  processing  within  a  Intel  Xeon  6-core  computer.  This  allowed  us  to 
obtain  a  six-fold  reduction  in  processing  time.  Eor  example,  the  time  required  for  the 
complete  SER  curve  over  the  same  range  of  INRs  was  reduced  from  3  h  to  0.5  h.  The  time 
reduction  allowed  us  to  perform  multiple  trials  within  a  more  realistic  length  of  time.  Eor 
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instance,  100  trials  could  be  completed  in  50  h  with  parallel  computing  instead  of  300  h 
without  parallel  computing.  Despite  the  six-fold  decrease  in  processing  time  achieved  with 
parallel  computing,  the  long  processing  time,  nevertheless,  imposed  a  practical  limit  on  the 
simulations. 

For  the  rest  of  this  chapter,  the  average  SERs  shown  are  calculated  over  100  trials.  The 
error  bars,  if  shown,  show  the  range  covered  from  the  5th  percentile  to  95th  percentile  SER 
values. 


4.2  Effect  of  Number  of  Coefficients  in  Channel  Estimate 

The  channel  estimator  in  the  preprocessor  provides  the  estimated  CIR  h{n)  through  an 
inverse-EET  operation.  The  number  of  coefficients  in  h(n)  is  determined  by  the  size  of  the 
EET.  In  most  practical  implementations,  not  all  the  coefficients  are  significant.  The  non¬ 
significant  coefficients  contribute  to  noise  [25].  In  Eai’s  implementation  of  the  channel 
estimator,  the  number  of  significant  coefficient  was  assumed  to  be  at  most  165,  and  the 
remaining  coefficients  were  regarded  as  noise  [11].  Since  the  channel  estimator  is  derived 
from  her  work,  the  channel  estimator  in  this  study  works  in  a  similar  fashion. 

We  first  investigate  the  effect  of  the  number  of  coefficients  in  the  CIR  has  on  the  resulting 
SER.  The  true  transmission  channel  has  17  coefficients.  To  remove  the  effect  of  the  trans¬ 
mission  channel  completely,  the  channel  estimator  must  have  at  least  the  same  number  of 
coefficients  as  the  transmission  channel,  that  is,  17  coefficients. 

We  keep  the  first  k  coefficients  that  come  out  of  the  channel  estimator  and  truncate  the 
rest.  Since  the  channel  estimator  produces  at  most  165  coefficients,  k  is  selected  such  that 
0  <k  <  165.  In  Eigure  4.3,  we  show  the  SER  as  a  function  of  the  INR  for  different  values 
of  k.  When  k  =  11,  the  channel  estimator  has  exactly  the  same  number  of  coefficients  as 
the  actual  transmission  channel.  We  denote  this  scenario  as  lx.  When  k  =  85,  the  channel 
estimator  has  five  times  more  coefficients  than  the  actual  transmission  channel.  We  denote 
this  scenario  as  5x.  Similarly,  the  scenario  whereby  k  =  165  is  denoted  as  9x,  rounded  to 
one  significant  figure.  Results  show  that  the  SER  degrades  when  the  length  of  the  estimated 
CIR  h{n)  is  higher  than  that  of  the  true  transmission  channel.  In  Eigure  4.4,  we  plot  the 
mean  SERs  with  their  error  bars  for  each  value  of  k.  The  variation  in  the  data  tends  to 
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increase  for  higher  values  of  k. 

This  behavior  is  likely  due  to  the  limitation  in  the  estimation  proeess  for  the  CIR.  A  typieal 
estimated  CIR  made  by  the  ehannel  estimator  is  shown  in  Figure  4.5.  There  are  only  a  small 
number  of  signifieant  eoeffieients.  The  non-signifieant  eoeffieients  are  non-zero  due  to  the 
limitation  in  the  estimation  proeess.  Note  that,  when  too  many  eoeffieients  are  used,  the 
noise  in  these  non-zero  eoeffieients  are  aeeumulated  and  results  in  signifieant  errors  in  the 
estimated  television  signal.  As  a  result,  the  adaptive  noise  eaneeler  is  not  able  to  reeover  the 
eommunieation  signal.  This  result  also  shows  that  a  good  ehannel  estimator  is  important 
in  this  seheme.  The  ehannel  estimator  must  not  only  be  able  to  aeeurately  estimate  the 
eoeffieients  itself  but  also  must  be  able  to  remove  the  non-signifieant  eoeffieients  due  to  the 
estimation  proeess  [25]. 


4.3  Effect  of  Error  in  Channel  Coefficients 

In  this  seetion,  we  first  assume  there  is  no  error  in  the  length  of  the  estimated  CIR  h{n). 
The  only  souree  of  error  in  h{n)  is  due  to  errors  in  the  eoeffieient  values.  We  run  the 
simulation  for  5%,  10%,  and  20%  errors  in  the  estimated  ehannel  eoeffieients  and  observe 
the  eorresponding  SERs. 
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The  plots  of  SER  against  different  INRs  are  shown  in  Figure  4.6.  The  5th  and  95th  per¬ 
centiles  are  shown  as  error  bars  in  Figure  4.7.  In  some  instances,  the  SERs  are  equal  to  zero. 
To  allow  us  to  plot  such  values  onto  the  logarithmic  axis,  they  are  represented  as  10^^  in 
the  plot.  In  general,  the  SER  increases  as  the  percentage  of  error  in  the  channel  coefficients 
increases.  There  are  occasions  where  the  SER  appears  to  increase  when  INR  increases, 
such  as  from  20  dB  to  25  dB  in  Figure  4.7b  and  4.7c,  respectively.  This  is  attributed  to  the 
large  outliers  in  the  associated  SER.  In  particular,  note  that  for  5%  and  10%  error  at  greater 
than  10  dB,  there  appears  to  be  no  error  bar.  Note  that,  in  these  cases,  fewer  than  5%  of  the 
SERs  values  are  nonzero.  That  is,  at  least  95%  of  the  SERs  are  equal  to  zero.  The  outliers 
beyond  the  95th  percentile  increased  the  mean  SER.  The  median  SERs  for  these  cases  are 
zero.  The  distribution  of  SER  values  for  each  INR  illustrates  this  skewed  behavior  and  is 
presented  in  Appendix  C. 

The  SERs  at  the  highest  percentage  error  of  20%,  however,  remain  below  10^^  for  INRs 
greater  than  0  dB.  This  shows  that,  with  an  appropriate  channel  estimation  process,  the 
approach  can  recover  the  weak  communication  signal  with  a  good  SER. 
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(a)  Plot  of  mean  SER  against  INR  with  42  000  symbols  transmitted  per  INR  over  100  trials.  The 
lower  and  upper  ends  of  the  error  bars  show  the  5th  and  95th  percentiles,  respectively. 


(b)  Plot  of  SER  against  INR  with  42 000 symbols  transmitted  per  INR  over  100  trials.  Each 
cross  shows  a  single  observation  of  the  SER  at  a  particular  INR.  SER  values  equal  to  zero  are 
represented  as  10“^  on  the  logarithmic  axis.  The  numbers  in  italics  just  above  the  .x-axis  indicate 
the  number  of  trials  with  SERs  equal  to  zero  at  a  given  INR  level. 

Figure  4.2.  Plot  of  SER  against  INR  with  42  000  symbols  transmitted  per  INR  over  100  trials. 
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Figure  4.3.  Plot  of  mean  SER  against  INR  with  different  number  of  coefficients  k  in  the  estimated 
CIR  averaged  over  100  trials.  The  number  in  parenthesis  represents  the  increase  in  the  number 
of  coefficients  used  in  the  estimated  CIR  as  compared  to  the  true  transmission  channel. 
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Figure  4.4.  Plot  of  mean  SER  against  INR  for  different  number  of  coefficients  k  in  the  estimated 
CIR  averaged  over  100  trials  with  error  bars.  The  error  bars  show  the  range  covered  from  the 
5th  percentile  to  95th  percentiles.  The  j-axis  range  in  all  plots  is  kept  constant  for  comparison, 
(a)  Mean  SER  for  k=\l  (lx),  85  (5x),  and  165  (9x).  (b)-(d)  Mean  SER  and  90%  confidence 
intervals  for  listed  values  of 
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Figure  4.5.  Typical  plot  of  estimated  CIR  h{n). 


Figure  4.6.  Plot  of  mean  SER  against  INR  for  different  errors  in  estimated  channel  coefficients 
averaged  over  100  trials.  SER  values  equal  to  zero  are  represented  as  10“^  in  the  logarithmic 
axis. 
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Mean  SER  Mean  SER 


Figure  4.7.  Plot  of  mean  SER  against  INR  for  different  errors  in  estimated  channel  coefficient 
values  averaged  over  100  trials  with  error  bars,  (a)  Mean  SER  for  5%,  10%,  and  20%  errors  in 
channel  coefficients,  (b)-(d)  Mean  SER  and  90%  confidence  intervals  for  listed  errors  in  channel 
coefficients. 
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CHAPTER  5: 
Conclusions 


In  this  chapter,  we  give  a  summary  of  the  results  obtained  and  suggest  recommendations 
for  future  work. 


5.1  Summary 

We  implemented  a  Simulink  model  to  examine  the  performance  of  an  adaptive  noise  can- 
celer  to  recover  a  weak  4-QAM  communication  embedded  in  a  strong  DTMB  television 
broadcast  signal.  We  successfully  integrated  the  DTMB  signal  generator  with  the  com¬ 
munication  signal  generator,  transmission  channel,  the  preprocessor  and  the  adaptive  noise 
canceler  to  form  a  working  Simulink  model.  Noting  the  long  processing  time,  we  imple¬ 
mented  parallel  computing  and  achieved  a  six-fold  reduction  in  the  processing  time.  Using 
the  resulting  model,  we  tested  the  adaptive  noise  canceler  under  different  scenarios.  Re¬ 
sults  showed  that  the  adaptive  noise  canceler  required  a  very  accurate  channel  estimation 
process  to  recover  the  weak  communication  signal.  With  the  correct  number  of  channel 
coefficients,  the  adaptive  noise  canceler  could  cope  with  at  least  up  to  20%  error  in  channel 
coefficients  to  produce  SERs  of  the  order  of  magnitude  10 

5.2  Recommendations  for  Future  Work 

The  following  recommendation  for  future  work  is  derived  from  some  questions  arising 
from  this  study. 

•  This  study  considered  the  ability  of  the  adaptive  noise  canceler  to  recover  a  4-QAM 
communication  signal  embedded  in  a  strong  DTMB  television  broadcast  signal.  Pos¬ 
sible  future  work  could  study  the  performance  of  the  adaptive  noise  canceler  with 
other  television  broadcast  signals  e.g.,  ATSC  or  DVB-T. 

•  The  channel  estimator  used  in  this  study  is  only  one  of  the  possible  channel  estimator 
approaches  available  for  DTMB.  For  example,  Zhang  proposed  an  iterative  channel 
estimator  scheme  [21].  A  possible  future  work  could  investigate  whether  the  adaptive 
noise  canceler  would  work  with  other  types  of  channel  estimator. 
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•  The  long  processing  time  noted  in  the  study  was  partially  mitigated  using  parallel 
computing  across  multiple  cores  within  a  single  machine.  Future  work  could  attempt 
to  further  reduce  the  processing  time  by  distributing  the  computational  workload 
across  a  computing  cluster  instead  of  a  single  machine  or  to  implement  the  model  in 
a  compiled  language  such  as  C. 
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APPENDIX  A: 

Simulink  Models 


In  this  appendix,  we  show  the  detailed  Simulink  models  used  in  this  study. 


Figure  A.l.  Overall  Simulink  model.  The  detailed  models  for  the  communication  signal  modulator,  DTMB  modulator,  multipath 
channel,  preprocessor  and  adaptive  noise  canceler  are  shown  in  Figures  A. 2  to  A. 5. 


lYansmit  Filter 

Figure  A. 2.  Simulink  model  for  communication  signal  modulator. 


Constant 


Figure  A. 3.  Simulink  model  for  DTMB  modulator. 


Channel  Dynamics 


Figure  A. 4.  Simulink  model  for  transmission  channel.  It  is  possible  to  switch  the  type  of  trans¬ 
mission  channel  using  the  selector  switch. 
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Figure  A. 5.  Simulink  model  for  preprocessor  and  adaptive  noise  canceler. 
mean  square  filter  and  recursive  least  square  filter  using  the  selector  switch, 
communication  signal  demodulators  are  shown  in  Figures  A. 6  to  A. 8. 


It  is  possible  to  switch  between  normalized  least 
The  detailed  models  for  DTMB  demodulator  and 


Gain 


Channel  Estimator 


(a)  Overall  view. 


Constant 


(b)  Channel  estimator. 

Figure  A. 6.  Simulink  model  for  the  DTMB  demodulator  and  the  channel  estimator  within  the 
DTMB  demodulator. 


Receive  Filter 


Figure  A. 7.  Simulink  model  for  the  first  communication  signal  demodulator. 


Receive  Filter 


Figure  A. 8.  Simulink  model  for  the  second  communication  signal  demodulator. 
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APPENDIX  B: 
MATLAB  Code 


In  this  appendix,  we  show  the  MATLAB  eode  that  have  been  developed  for  this  study. 

B.l  Code  to  Run  Simulation 

runsim.m  is  a  funetion  that  initializes  the  environment  for  the  simulation,  runs  it  and  re¬ 
turns  the  results.  It  ealls  adapt_dtmb_mqam_init  .m  to  initialize  the  parameters  that  do  not 
vary  between  simulation  runs  and  initialize  the  parameters  that  do  vary  between  simulation 
runs. 


function  [ SER_totalerr ,  timetaken,  channelhl,  simOut]  =  ... 

runsim ( INR_dB,  stopsim_maxnumsymbols,  channelcoef f_error , 
selectf irstKcoef f s,  numTrials,  isSilent) 

%RUNSIM  Run  thesis  simulation. 

%  [ SER_totalerr ,  timetaken,  channelhl,  simOut]  =  ... 

%  runsim ( INR_dB,  stopsim_maxnumsymbols ,  channelcoef f_error, 

%  selectf irstKcoef fs,  numTrials,  isSilent) 

%  Returns 

%  simOut  Contains  results  of  simulation 

%  timetaken  Time  taken  for  this  run 


Inputs 

lNR_dB  Interf erence-to-Noise  ratio,  defined  as  comms  signal  over 

channel  noise. 

stopsim_maxnumsymbols  Run  simulations  for  specified  number  of 
symbols . 

channelcoef f_error  Error  in  estimation  of  channel  coefficients 
E.g  Specify  0.2  for  a  20%  error  in  coefficient, 
selectf irstKcoef fs  The  channel  estimator  attempts  to  model  the 
channel  as  a  FIR  filter  and  outputs  255  coefficients.  This 
parameter  selects  the  first  selectf irstKcoef fs  coefficients  to  be 
used  in  the  signal  subtraction. 
numTrials  Number  of  trials  to  repeat 

isSilent  If  true,  the  Simulink  model  is  loaded  silently. 
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%  No  window  will  be  seen.  This  option  is  used  most  commonly  when 
%  doing  batch  runs. 

o, 

o 

%  See  also  MAINPAR  and  SIM. 


o,  o_ 
o  o 

f ilenameModel  =  ' adapt_dtmb_mqam ' ; 
filenamelnit  =  strcat (f ilenameModel,  '_init'); 

if  isSilent 

%  Load  Simulink  model  w/o  opening  a  window 
load_system (filenameModel) 
else 

%  Open  Simulink  model 
open_system (filenameModel) 
end 

%%  Initialize  the  model  variables 
%  Run  initialization  script 
eval (filenamelnit) 

%  Set  the  power  of  noise  generator  based  on  INR_dB 
sigma_w  =  10'' (  (-SIR_dB  -  INR_dB)  /  10); 

%  Set  the  value  of  selectf irstKcoef f s 

channel . selectfirstKcoeffs  =  selectfirstKcoef f s; 

%  if  isRapid 

%  paramNameValStruct . SimulationMode  =  'rapid'; 

%  end 

%  Set  options  for  running  sim  command 

%  See  MATLAB  help  for  more  details  on  parameters  for  sim  command 
paramNameValStruct . ReturnWorkspaceOutputs  =  'on'; 
paramNameValStruct . SrcWorkspace  =  'current'; 

SER_totalerr  =  zeros (  numTrials,  3  ); 

channelhl  =  zeros (  numTrials,  length (channel . hi ) *1  ) ; 
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o,  g_ 
0  0 


Run  simulation 


timestart  =  tic; 

for  iTrials  =  1  :  numTrials 

%  Re-run  initialization  script  to  randomize  the  txn  channel 
eval { f ilenamelnit ) 

simOut  =  sim (filenameModel,  paramNameValStruct ) ; 
SER_totalerr (iTrials,  :)  =  simOut . get { ' SER_totalerr ' ) ; 
channelhl ( iTrials ,  :)  =  channel. hi  ; 

end 

timetaken  =  toe (timestart ) ; 


%%  Compute  results 

%  signalcomms_power  =  mean (abs (reshape (.. . 

%  squeeze (simOut . get ( ' signalcomms ' ) ) , 1 , [ ] ) ) • ^2 ) ; 

%  signaltv_power  =  mean (abs (reshape (.. . 

%  squeeze (simOut . get ( ' signal tv ')),!,  []))  .^2); 

%  signalnoise_power  =  mean (abs (reshape (.. . 

%  squeeze (simOut . get ( ' signalnoise ' )  ) , 1 ,  [  ] ) )  •  ^2 ) ; 

%  Returns  SER,  total  errors,  symbols  txn 
%SER_totalerr  =  simOut . get (' SER_totalerr ') ; 

%  Prints  out  signal-to-interf erence  ratio  (SIR)  and 
%  interf erence-to-noise  ratio  (INR) 

%fprintf('TV  to  Comms  Ratio  =  %0.2f\n',  ... 

%  10*logl0 (signaltv_power/signalcomms_power)  ) 
%fprintf (' Comms  to  Channel  Noise  Ratio  =  %0.2f\n',  ... 

%  10*logl0 (signalcomms_power/signalnoise_power)  ) 


%  Prints  out  time  taken 

fprintf ( ' Time  taken  =  %0.2f  minin',  timetaken/ 60 ) 
%  Prints  out  SER 

%fprintf ( ' SER  =  %0.2f  \n',  SER_totalerr ( 1 )  ) 

fprintf ( ' \n ' ) 
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end 


B.2  Code  for  Initialization 

adapt  _dtmb_mqain_init .  m  is  an  helper  funetion  that  contains  the  initialization  parameters 
for  the  simulation.  Parameters  that  are  meant  to  be  varied  between  simulation  runs  are 
initialized  in  runsim .  m. 

%  ADAPT_DTMB_MQAM_INIT  A  helper  function  to  put  the  initialization  of 
%  parameters  in  a  single  location.  It  is  not  meant  to  be  executed 
%  directly  by  the  user. 

%%  DTMB  code  starts  here 

%  DTMB  Parameters  for  different  cases: 
dtmb.M  =16;  %  M-QAM 
dtmb.EbNO_dB  =  15;  %in  dB 

%  Rayleight  Multipath  CH  with  Doppler: 

%  doppler  freq.  in  Hz  (change  when  testing  under  different  environment) 
dtmb . f D  =  80; 

%  Multipath  time  delay  in  (sec) 

dtmb . DelayPath  =  [0  0.2e-6  0.5e-6  1.6e-6  2.3e-6  5e-6] ; 

%  Multipath  power  in  (dB) 
dtmb.Pwd  =  [-3  0  -2  -6  -8  -10]; 

%  Frequencies 

Fs  =  7.56e6;  %  sampling  and  symbol  rate  in  Hz 

dtmb.NFFT  =  3780;  %  FFT  length 

%  compute  noise  variance  according  the  specified  Eb/NO 
dtmb.EbNO  =  10^ (dtmb . EbN0_dB  /lO); 

%  Received  delay  setting  in  error  rate  calculation  block 
dtmb . TD  =  zeros  (1,16); 

%  M=16,  16-QAM 
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dtinb.TD{16)  =  30244; 

%  M=4,  QPSK 
dtinb.TD{4)  =  15122; 

%  M=2,  BPSK 
dtinb.TD{2)  =  75  60; 

%  noise  variance 

%  signal  power=unity  power=OFDM  data  power 
dtmb.Ps  =  1; 

%  prefix  length  (420) 
dtmb . L  =  dtmb . NFFT/ 9 ; 

%  OFDM  data  length (3780) 
dtmb.N  =  dtmb. NFFT; 

%  prefix  power  boost  factor 
dtmb . r  =  2 ; 
dtmb.eps  =  eps; 

%  PN  sequence  for  GI 
%  QPSK  signal 

dtmb.pO  =  sign (randn (1, 255) )  +  li*sign (randn (1, 255)  )  ; 

%  normalized  prefix  power  to  twice  of  the  OFDM  data  part 
dtmb.pO  =  dtmb .pO/ (sqrt (dtmb .NFFT) )  ; 

dtmb.pre  =  dtmb .pO (255-83+1 : 255) ; 

dtmb. post  =  dtmb . pO ( 1 : 82 )  ; 

dtmb.pn  =  [dtmb.pre,  dtmb.pO,  dtmb. post]; 

dtmb.NGI  =  255+82+83;  %  Length  Guard  Interval 

dtmb.PnO  =  f f t (dtmb . pn ( 1 66 : dtmb . NGl ) )  ; 

dtmb.PnOinv  =  con j (dtmb . PnO ) . / (abs (dtmb . PnO ) . ^2  +  dtmb.eps); 
dtmb.zeropn  =  [ zeros ( 1 ,  dtmb . NFFT)  ,  dtmb.pn]; 

%  bit  rate 

dtmb.Fb  =  log2 (dtmb . M) * (dtmb . NFFT/ (dtmb . NFFT+dtmb . NGI )) *Fs ; 


0,0,  g,  g,  g, 

0  0  o  o  o 

%  Code  added  to  make  DTMB  model  work  in  comm  signal 
%  This  parameter  is  set  in  the  runsim.m  file 
%sigma_w=0 . Ole-6;  %  Noise  Variance 
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%%  Dynamics  of  transmission  channel 
%  with  multipath 

channel. PdB  =  [0,  -10,  -15];  %  Power  Profile  in  dB ' s 

channel. tau  =  [0,  10/Fs,  15/Fs];  %  time  delays  in  sec 
channel. n  =  round { channel . tau*Fs )  ;  %  time  delays  in  samples 

%  convert  indices  from  zero-based  to  one-based 
channel. n  =  channel. n  +  1; 

channel. hO  =  (randn (size (channel . PdB) )  +  ... 

li*randn (size (channel .PdB) ) )  /  sqrt(2); 

channel . sqrtP  =  10 . ^ (channel . PdB/20 ) ; 
channel. hi  =  zeros(l,  max (channel . n) )  ; 
channel . hi (channel . n)  =  channel. hO  .*  channel . sqrtP ; 

%  normalize  for  unit  power 

channel. hi  =  channel. hi  /  sqrt ( (channel .hl*channel .hi ')) ; 
channel. hi  =  [ 0 , channel . hi ] ;  %  strictly  causal 


%  Estimated  Channel 

%  with  random  error  in  coefficients  as  specified 
%  by  channelcoef f_error 

channel . hlest  =  channel. hi  +  channelcoef f_error  *  channel. hi  .*  ... 

(  randn ( size (channel . hi ) )  +  li*randn (size (channel .hi) )  )  /  sqrt(2) 

%  channel . hlest  =  channel. hi;  %  assume  no  errors  first 

%  This  parameter  is  set  in  the  runsim.m  file 
%  channel . selectfirstKcoeffs  =  17; 

%%  Paramters  for  the  Weak  Comms  Signal 
comms.Fc  =  Fs/10;  %  carrier  frequency 

comms. Fdw  =  comms. Fc/10;  %  data  rate 

%  comms. M  Number  of  symbols  in  comms  signal  i.e.  4-QAM  or  16-QAM. 
comms. M  =4;  %  4-QAM  comms  signal 

0,0,  g,  o,  o, 
o  o  o  o  0 
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%  Signal  to  Interference  Ratio  (SIR) 

%  defined  as  TV  signal  /  comms  signal 
SIR_dB  =  30; 

%  LMS  weights  vector 

%n  =  0 : length (dtmb . DelayPath) -1 ; 

LMSf liter . length  =  50; 

LMSfilter .weights  =  zeros ( 1 , LMSf liter . length) ; 
LMSf liter .mu  =  0.0001;  %  stepsize  of  LMS 

%  Number  of  data  points  to  keep  for  debug 
keepLastMpoints  =  1; 

9-  9-  2-  9- 

o  o  o  o 


B.3  Code  to  Enable  Parallel  Computing 

mainpar .  m  is  a  script  that  distributes  the  running  of  stimulation  runs  aeross  multiple  pro- 
eessors,  allowing  parallelization  to  take  plaee. 

%  MAINPAR  Main  job  to  start  other  tasks  for  parallerization . 

%  This  is  the  main  file  to  run  to  generate  a  complete 
%  SER-INR  curve.  To  determine  the  SER  for  a  single  INR  value, 

%  see  RUNSIM. 

o, 

o 

%  The  following  internal  variables  within  the  script  file 
%  control  the  parameters  of  the  simulation. 

o, 

o 

%  INR_dB  is  an  array  that  sets  the  INRs  to  run  for. 

%  stopsim_maxsym  sets  the  number  of  symbols  transmitted 
%  for  that  particular  value  of  INR  before  stopping  the 

%  simulation. 

%  num_Trials  is  the  number  of  trials  to  repeat. 

%  selectf irstKcoef f s  sets  the  number  of  coefficients  to  keep 

%  in  the  estimated  channel  impulse  response. 

%  channelcoef f_error  sets  the  percentage  error  in  the  channel 

%  coefficient  values. 


37 


%  Note  that  the  Simulink  model,  adapt_dtmb_mqam . six,  determines 

%  whether  the  simulation  is  using  selectf irstKcoef f s  or 

%  channelcoef f_error . 

o, 

o 

%  See  also  RUNSIM,  PARFOR  and  MATLABPOOL. 

%clearvars 
%close  all 

%%  Set  parameters  here 
INR_dB  =  (-5  :  5  :  50)  '  ; 

%  start  with  small  number  to  test  parfor 
stopsim_maxsym  =  4.2e3  *  10  *  ones (size {INR_dB) )  ; 

%  Adjust  channelcoef f_error  or  selectf irstKcoef f  depending  on  test  case 

%  Remember  to  adjust  the  Simulink  model  and  the  filename  to  save  the 

%  results  to  as  well 

channelcoef f_error  =  0.3; 

selectf irstKcoef fs  =  17; 

numTrials  =  100; 

%  length  of  channelhl .  Depends  on  Fs  and  channel  parameters 
%  this  is  used  to  do  preallocation 
channeltxn_len  =  17; 

%%  Loop  thru  simulations 


elapsed  =  zeros ( size ( INR_dB) ) ; 

SER  =  zeros (  numTrials,  3,  size(INR_dB,  1)  ); 

channelcoef fs  =  zeros (  numTrials,  channeltxn_len*l ,  size(INR_dB,  1)  ); 

parfor  1=1:  length ( INR_dB) 

%[SER{i,:),  elapsed(i),  ~]  =  runsim (INR_dB (1) ,  stopsim_maxsym ( 1 ) , 

%  comms_M,  selectf irstKcoef fs ,  numTrials,  true); 

[SER{ elapsed(i),  channelcoef fs i)  ,  ~]  =  ... 
runsim ( INR_dB ( 1 ) ,  stopsim_maxsym { i ) ,  ... 

channelcoef f_error,  selectf irstKcoef fs ,  numTrials,  true)  ; 
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end 


%%  Tabulate  results 

%results  =  [INR_dB  SER]  ; 
results  =  [  ... 

reshape!  repmat {INR_dB, 1, numTrials)  ' ,  [numTrials*length ( INR_dB)  ,  1]), 
reshape!  permute ! SER, [132])  ,  [numTrials*length !INR_dB)  3]  ) 

]  ; 


channeltxn  =  reshape!  permute !channelcoeffs,  [13  2])  ,  ... 

[  numTrials*length !INR_dB) ,  channeltxn_len]  ) ; 

%  Save  handle  of  figure 
hdl  =  zeros ! 1 , 2 )  ; 
hdl  !  1 )  =  figure ; 

plot  !  results  !:, 1 ) ,  results]:, 2)  ) 

xlabel ! ' INR  !dB) ' ) 

ylabel !  sprintf !  'SER  over  %d  trials',  numTrials)  ) 

title!  sprintf  !' selectfirstKcoeffs  =  %d  channelcoeff error  =  %.2f',  ... 

selectfirstKcoef fs,  channelcoef f_error )  ) 

results_cumu  =  [  INR_dB,  zeros!  length ! INR_dB) ,  4 ) ] ; 

for  i  =  1 : length !INR_dB)  , 

%  Sum  of  error  encountered  and  symbols  transmitted 
results_cumu !i,  4:5)  =  ... 

sum!  results!  numTrials* !i-l) +1 :numTrials*i,  3:4  )  )  ; 

%  Mean  SER 

results_cumu [i,  2)  =  ... 

mean!  results!  numTrials* !i-l) +1 :numTrials*i,  2  )  ); 

%  Median  SER 
results_cumu [i,  3)  =  ... 

median!  results!  numTrials* !i-l) +1 :numTrials*i,  2  )  ); 

end 

hdl [2)  =  figure; 

plot ! results_cumu ! : , 1 ) ,  results_cumu ! : , 2 ) ,  ... 

results_cumu [:,!),  results_cumu ! : , 3 )  ) 
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xlabel { ' INR  (dB) ' ) 

ylabel {  sprintf  (  ' SER  over  %d  trials',  numTrials)  ) 

title (  sprintf (  strcat (' Cumulative .  selectf irstKcoef f s  =  ',  ... 

' %d  channelcoeff error  =  %.2f'),  ... 

selectfirstKcoef fs,  channelcoef f_error )  ) 

legend {' Mean ' ,  'Median') 

hdl ( 3 )  =  figure ; 

semilogy {results_cumu ( : , 1 ) ,  results_cumu { : , 2 ) ,  ... 

results_cumu {:,!),  results_cumu ( : , 3 )  ) 

xlabel { ' INR  (dB) ' ) 

ylabel {  sprintf (  'SER  over  %d  trials',  numTrials)  ) 

title (  sprintf (  strcat (  'Cumulative,  selectf irstKcoef fs  =  ',  ... 

' %d  channelcoef f error  =  %.2f'),  ... 

selectfirstKcoef fs,  channelcoef f_error)  ) 
legend {' Mean ' ,  'Median') 

%%  Save  results 
f ilename_results  =  ... 

sprintf (  strcat ( ' Results/ results_%s_NLMS_selectfirstKcoef f%d_' 
' %dtrials_%dsymbols ' )  ,  ... 

datestr (now,  ' yyyymmdd ' ) ,  selectfirstKcoef fs,  numTrials,  ... 
stopsim_maxsym ( 1 )  ); 

save_results ( f ilename_results ,  . . . 

hdl,  results,  results_cumu,  channeltxn,  ... 
numTrials,  selectf irstKcoef fs,  channelcoeff_error) ; 

%%  Clean  up  variables 

%clear  INR_dB  stopsim_maxsym  selectf irstKcoef fs  numTrials  1  ... 

%  channelcoef f_error  channelcoef fs  channeltxn_len  SER 

clear  INR_dB  stopsim_maxsym  1  channelcoef fs  channeltxn_len  SER  hdl 

fprintf('All  done!\n') 


function  save_results (f ilename_main,  handle,  results,  ... 

results_cumu,  channeltxn,  numTrials,  selectfirstKcoef fs, 
channelcoef f_error  ) 

%  SAVE_RESULTS  Saves  the  results  of  a  completed  simulation  run 

o 

0 
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%  save_results  ( f  ilenaine_main,  handle,  results,  ... 

%  result s_cumu,  channeltxn,  numTrials,  selectf IrstKcoef f s ,  . 

%  channelcoef f_error  ) 

o, 

o 

%  This  function  should  not  be  called  directly  by  the  user  under 

%  normal  situations.  Designed  to  be  a  helper  function  to  MAINPAR 

o, 

o 

%  See  also  MAINPAR. 

%  Saves  the  variables 

save ( f ilename_main,  'results',  ' results_cumu '  ,  'channeltxn',  ... 
'numTrials',  ' selectf IrstKcoef fs ' ,  ' channelcoef f_error '  ) 

%  saves  Figure  1,  2  and  3 

%  Intentionally  not  use  '!'  as  loop  counter  to  avoid  clash 
for  iCount  =  1:3 

filename  =  strcat  (  f ilename_main,  sprintf  (  '-%d',  iCount  )  ); 

saveas (handle ( iCount )  ,  filename,  'fig') 
end 

%  Clean  up  variables 

clear  f ilename_main  filename  iCount 
end 
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APPENDIX  C: 
Supplementary  Plots 


In  this  appendix,  we  present  additional  plots  to  supplement  the  information  provided  in 
Chapter  4. 
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Figure  C.l.  Plot  showing  distribution  of  the  100  observations  of  SER  for  each  INR  value  for 
different  number  of  coefficients  k  in  h{n).  The  y-axis  range  is  kept  constant  throughout  for 
comparison. 
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Figure  C.2.  Plot  showing  distribution  of  the  100  observations  of  SER  for  each  INR  value  for 
different  percentages  of  error  in  the  channel  coefficients.  The  j-axis  range  is  kept  constant 
throughout  for  comparison.  SER  values  equal  to  zero  are  represented  as  10“^  on  the  logarithmic 
axis.  The  numbers  in  italics  just  above  the  .x-axis  denote  the  number  of  SER  values  that  are 
equal  to  zero. 
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